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ABSTRACT 



A sample of 10 nearby intermediate-type active galactic nuclei (AGN) drawn from 
the Sloan Digital Sky Survey (SDSS-DR7) is presented. The aim of this work is to pro- 
vide estimations of the black hole mass for the sample galaxies from the dynamics of 
the broad line region. For this purpose, a detailed spectroscopic analysis of the objects 
was done. Using BPT diagnostic diagrams we have carefully classified the objects 
as true intermediate-type AGN and found that SO%'^ljf^^ are composite AGN. The 
black hole mass estimated for the sample is within 6.54+0.16 < log Mbh < 7.81+0. 14. 
Profile analysis show that five objects (J120655.63+501737.1, J121607.08+504930.0, 
J141238.14+391836.5, J143031. 18+524225.8 and J162952.88+242638.3) have nar- 
row double-peaked emission lines in both the red (Ha, [N n]/lA6548,6583 and 
[S II]/l/16716,6731) and the blue (Bfi and [O III]^/14959,5007) region of the spec- 
tra, with velocity diff"erences (AV) between the double peaks within 1 14 < AV < 256 
km s-\ Two of them, J121607.08+504930.0 and J141238.14+391836.5 are candi- 
dates for dual AGN since their double-peaked emission lines are dominated by AGN 
activity. In searches of dual AGN; Type 1, Type II and intermediate-type AGN should 
be carefully separated, due to the high serendipitous number of narrow double-peaked 
sources (50%+ 14.4%) found in our sample. 

Subject headings: galaxies: active - galaxies: Seyfert -galaxies: nuclei-quasars:emission lines 
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1. Introduction 



Among nearby active galactic nuclei (AGN) the so-called Seyfert ( Sy) galaxies are common 



found . Seyfert galaxies properties were described in the early work by ^Khachikian & Weedman 



19741) who classified them as Sy 



optical spectra. A few years later, 



and Sy 2 depend ing on the line widths observed in their 



Osterbrockl (|1981|) introduced different sub-classes of Sy 



galaxies according to the relative intensity showed by the broad emission lines in their spectra. In 
particular, he introduced the intermediate Sy class, spanning from Sy 1.2 to Sy 1.9. For instance, 
a Sy 1.8 galaxy shows weak, but readily visi ble broad Ha and H/j emission lines, while a Sy 1.9 



shows only a broad Ha component. Recently, 



Tadhuntei 



(|2008b reviewed the general classification 



and unification of AGN where Sy 1 and Sy 2 galaxies are more generally grouped as Type 1 and 
Type 2 AGN. Since intermediate-type Sy galaxies do show broad emission line components, they 
belong to the Type 1 class, but they are more diflicult to find since sometimes the intensity of 
these components can be very low. We will refer to these objects as intermediate-type AGN to 
distinguish them from the more general Type 1 class. 

Early studies on intermediate-type AGN, have proposed that the origin of spectroscopically 



observed weak emission lines pr oduced in the clouds of the Broad 



explained by internal reddening (lOsterbrock 



based in the unified model for AGN dUrry & Padovani 



1981 



Goodrich 



1990 



ine R egion (BLR) could be 



19951). Another explanation 



19951) assumes that a dusty-torus located 



some parsecs away from the central source blocks part of the emission. Other possibility is that 
their weak lines could be due to an intrinsic ally faint non-therm al continuum emission, maybe due 



to variations in the ionizing continuum (see 



Trippe et al. 



2010l and references therein). 



On the other h and, the availabi 



Sky Survey (SDSS; 



York et al. 



ity of high-quality spectra of AGN in the Sloan Digital 



20001) makes it possible to identify intermediate-type AGN and 



furthermore study th em as a c 



for nearby AGN (see 



Ho 



ass, in order to find their place in the black hole demography 



2004|) . In general, new estimates of fundamental parameters like the 
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black hole (BH) mass or Mbh in AGN can contribute to the understanding of the associated 
fueling mechanisms. The M^h can be estimated using different empir ically calibrated sca ling 



relations in samples designed to study particular types of AGN (e.g. 



Wandel et al 



Vestergaard & Peterson 



19991 : 



2006|) . In particular, BH ma ss estimates can be obta ined using the 

2 



200( 



2005h and BLR radius. 



correlation between the AGN continuum luminosity (|Kaspi et al. 
By combining continuum luminosity with the width of the Hfi emission line (in objects that 
have z<l, otherwise this line goes to the IR) it is possible to estimate the mass of the black 
hole by using a single AGN s pectrum, yielding BH masses that are accurate to a factor of ~3 



(|Vestergaard & Peterson 



2006h 



In the case of intermediate-type AGN, estimates of the Mbh need to be done systematica" 



Previous works in this type of AGN are ba s ed on single or very few objects (e.g 



1981 



Goodrich 



Ouillen et a 



Trippe et al 



1989; 



2000 



Rafanelli et al 



Torrealba et al 



1993 



Xanthopoulos 



19961 : 



Osterbrock 



Popovic & Mediavilla 



1997 



ly. 



20061) . A larger amount of objects was recently studied by 



(|2010f) . These authors selected a sample of 34 previously classified Sy 1.8/1. 9 and 
found that 18 are missclassified objects. Therefore, it is necessary to carefully classify candidates 
in order to establish if they conform a true sample of intermediate-type AGN. 

Among the Type 1 and Type 2 AGN samples drawn from the SDSS, the sample of 
double-peaked narrow emission line AGN stands out. Usually they show [O III] line profiles 
with velocities splitting in the range 151 to 1314 kms~^ The relevance of studying these kind 
of AGN is based on th e idea that they are candidates for sub-kpc or kpc-scale dual AGN (e.g. 



Comerford et al 



20091). Discovering new dual (or binary, if their nuclei appears separated on 



pc-scales) AGN candidates is fundamental to constrain models of galaxy formation and evolution. 
Since it is commonly believed that mergers can trigger or enhance nuclear accretion, it is 
fundamental to establish their frequency by means of detecting dual or binary AGN. Nevertheless, 
it is a fact that direct observational evidence for binary super massive black holes still remains 
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scarce (see lMcGurk et al.ll201 iL and references therein) 



Double-peaked [O in 



2009 



Rosario et al. 



20101 : 



line profiles could arise from sev eral mechanisms (e.g. JXu & Komossa 



Fu et al 



2011al : 



Shen et al. 



201lh : (1) the orbital motion of a binary 



AGN, (2) gas kinematics in extended narrow-line regions, and ( 3) unreso l ved nu clear gas 



Fu et al. 



( 2011b ) have obtained 



kinematics (e.g., aligned outflows or disk rotation on small scales), 
high-resolution imaging of 106 narrow double-peaked AGN (37 Type 1 and 69 Type II) and found 
that 31 have companions within 3 ". Their integral field spectroscopic observations has allowed 
them to group the narrow double-peaked AGN in the above three broad categories, according 
to the origin of the double-peaked [OIII] line profiles, finding that only 4.5-12 % are binaries. 



Other studies on the fre quency of these objects among AGN (e.g.. 



2010 



Smith et al. 



Wang et al. 



2009 



Liu et al 



2010h have show that there are 340 objects between 0.008 < z <0.686 or. 



equivalently, ~1 % of t he entire SPSS AGN sample harboring possible dual AGN. Nevertheless, 



a recent study done by 



Fu et al. 



(|2012|) shows that several objects classified in the literature as 
"binary" are, most likely, single AGN with extended narrow-line regions. Therefore, the frequency 
of double-peak objects and binaries among intermediate-type AGN are still unknown, and so this 
is worth to investigate. 

The aim of this work is to isolate intermediate-type AGN objects and estimate their black hole 
mass. In particular, we are interested in providing new estimates that can contribute to the study 
of the black hole demography in nearby AGN. For this purpose, we present the spectral analysis 
of a sample of 10 nearby (10'^^ < Lboi < lO'^'^erg s"') intermediate-type AGN (i.e. intermediate 
Sy galaxies) . Spectral analysis will provide a more precise classification through the use of the 
Baldwin-Phillips-Terlevich (BPT) diagnostic diagrams. A detailed profile decomposition allowed 
us to discover new narrow double-peaked AGN candidates among the objects in the sample. 

The paper is organized as follows: In Sect. |2]we describe the sample selection and in Sect. [3] 
the analysis of the SDSS spectroscopic data. Sect. |4] shows the results obtained for the sample and 
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relevant properties of some individual objects. The discussion and conclusions are presented in 
Sect. [51 An Appendix on the available X-ray data for one of the selected sources is presented. The 
cosmology adopted for this work is Hq =70 km s"^ Mpc"\ Q.^ =0.3 and O.^ =0.7. 



2. Sample Selection 

We selected objects flagged as quasars (QSO) in the SDSS database that have Petrosian 
g-band magnitudes between 14 < < 17 and redshifts within 0.026 < z < 0.12, i.e., a sample 
of nearby AGN. From this sample, we chose only those galaxies that have spectroscopic data 
in the SDSS database in order to accurately obtain their spectral classification. An additional 
condition imposed was that visually all selected galaxies should have spectral characteristics of 
intermediate-type AGN. We ended up with a sample of 23 candidates that were photometrically 
observed with the Nordic Optical Telescope (NOT). Results obtained with these data are presented 
in a companion paper (Benitez et al. 2012b, hereafter Paper II). A final sample of 10 candidates for 
intermediate-type AGN were chosen for this work, these are presented in Table [1] The remaining 
objects were taken out from the sample since after a preli minary analysis of their spe ctra, we 



found that some of them were either narrow line Seyfert 1 dOsterbrock & Pogge 
1 objects. 



1985b or Seyfert 



3. Analysis of SDSS-DR7 spectra 



The optical SDSS-DR7 spectra of our sample were processed using the STARLIGHT 



code (|Cid Fernandes et al. 



2005 



Mateus et al- 



to study some AG^ 



Cid Fernandes et al. 



samp 



2011 



es (see 



Bian et al. 



Coziol et al 



2011 



2007 



20061). This code has been successfully used 



Stasihska et al. 



Mezcua et al 



20111 : 



2006; 



Rafanelli et al. 



Leon-Tavares et al. 



first run, all spectra were corrected for Galactic extinction using the dust maps of 



20091 : 



2011 



)■ Ina 



Schlegel et al. 
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(Il998|) and the extinction law of ICardelli et alj (Il989|) . Also in this first run, in order to 
obtain the signal-to-noise ratio (S/N) we have masked the most prominent ernission lines that 



are known to be present in the composite spectra of QSOs (see 



Francis et al. 



1991 ). So, the 



following lines were masked: [O II]i^3726,3729, [Ne III]/{3869, H6i3970, H5/14102, Hy A4340, 
li/3A4S6l, [O III]i/l4959,5007, He U5876, [O I]/16300, HQr/16562, [N II]/l/J6548,6583, and 
[S II]/ii67 17,6731; together with the ISM absorption from NaDi5890 A). The signal-to-noise 
ratio (S/N) was obtained for each spectrum in the wavelength range of 4730 to 4780 A and is 
presented in Col. 2 of Table [2l 

Based on the S/N previously obtained, each spectrum was randomly perturbed 100 times 
using a Monte Carlo (MC) seed. Our lowest S/N value was 5.9, so the perturbations varied +16% 
of the signal value. We obtained 100 simulated spectra per object that were also analyzed with 
STARLIGHT, masking again the prominent emission lines as we did in the first run. From these 
simulated spectra, we obtained a collection of 100 models for the absorption spectrum of the host 
galaxy, from which we will obtain an estimation of the velocity dispersion. We chose the median 
of the distribution of models as our best model of the host-galaxy absorption spectrum. The stellar 
velocity dispersion cr^ was estimated from the median value of the distribution and the quoted 
errors are the rms values (see Table [5]). In all cases the obtained o^^ , is larger than 65 km s \ 



Liu et al. 



20101). The AGN continuum flux 



which is the nominal SDSS velocity resolution (e.g. 
at 5100 A was estimated from the the collection of 100 power-law fittings. Again we chose the 
best fit from the median value of the distribution of models to obtain F 5100 A and the rms as 
the associated error. This flux will be used to derive the BH mass, as is described in Sect. 14.11 
Therefore, STARLIGHT provided us with the best synthetic spectra for our AGN sample. 

The original SDSS spectra and the final pure AGN spectra, together with the obtained fits 
for the sample objects are shown in Figs.[I]to[5l These figures show how the stellar contribution 
is subtracted with STARLIGHT from the composite spectra. This is essential for our study since 
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the absorption synthetic spectra were used to estimate the velocity dispersion and the pure AGN 
spectra to perform a line profile analysis. An example of the best model obtained for object #9 is 
shown in Figure [6l where it can be seen how STARLIGHT is capable to fit the absorption lines 
used to estimate the velocity dispersion along with the stellar continuum. The absorption lines 
used to estimate cr* are marked in the blow-up spectrum. 

The complete optical spectroscopic parameters obtained for the sample are presented in 
Table 121 The FWHM for both the Ha and HyS broad components have errors <10%. The line ratios 
have also errors <10% except for object #10 that has errors ~18%. The profile decomposition in 
all cases was done using the algorithm PeakFitQ. From our profile decomposition analysis, we 
find that in five objects (#1, #3, #4, #5 and #8) the [O III]/l^4959,5007, [S II]i.l67 16,6731, [N II] 
/1/16548,6583 doublets and the narrow components of Ha and HJ3 clearly show a double-peaked 
profile. The spectral profile decomposition of these five narrow-double peaked AGN are shown in 
Figs. |7]to[lT] We have checked the redshifts using some absorption lines from the host galaxy 
spectra and found that they are in good agreement with the values reported in NEeP, which are 
based mainly on emission lines. Therefore, the Gaussians centroids were forced to correspond to 
the host galaxy redshift. For all objects, the widths of the Gaussians fitted to the narrow lines have 
the same Doppler broadening, and the intensity ratio [O III]/l5007/[O III]/14959 was fixed to its 
theoretical value. Table [3] presents the velocity difference defined as AV = Vred - Hiue for each 
object with double-peaked narrow emission lines for various lines. 



Baldwin et al, 



(|198l[ ) propose a suite of diagnostic diagrams (known as BPT) to classify 



emission-line galaxies in order to disentangle the dominant energy source. The diagrams are based 



'Systat Software Inc. http : //ww w . sigmaplot . coin| 

^The NASA/IPAC Extragalactic Database (NED) is operated by the let Propulsion Labora- 
tory, California Institute of Technology, under contract with the National Aeronautics and Space 
Administration. 
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on fou r optical line ratios [O III]/H/3, [N llj/Ua, [S II]/Ha and [O I]/Ha. Following iKewlev et al 



dlOOq) . we present in Figure [T2ltwo BPT diagnostic diagrams that show the location of the sample 



objects. These diagrams used line ratios [O III]/HjS, [N II]/Hq' and [S II]/Ha. 

We note that for the narrow double-peaked objects we present the obtained line ratios 
separately using the blue and red components. This allowed us to see how objects can change 
their spectral classification in the diagnostic diagram. We have found that 8 objects are composite 
AGN (Starburst plus AGN). In Table [5] we present the spectral classification obtained for all 
objects. In some objects we obtained different classifications since they appear in different places 
in the two BPT diagrams. In particular, object #2 appears as a Starburst (SB) galaxy in one of 



the BP' 



Whittle 



diagrams. In addition, since all objects present a broad line component, and following 



(119921 ). we have obtained their Sy type based on Fsoot/F//^ = R. Then, we have assigned 
Sy 1 (R<0.3), Sy 2(0.3 <R<l),Sy 1.5(1<R<4), Sy 1.8 (R> 4), and Sy 1.9 if only broad Ha 
was seen. 



4. Results 



4.1. Mbh estimates 



The Mbh for all obj ects in the sample were indirectly derived using an empirically calib rated 



photoionization method ( Wandel et al. 



19991 : 



Vestergaard.2002i : iVestergaard & Peterson 



derive d BH masses using the SDSS-DR7 spectra using equation 5 from 



2006). We 



Vestergaard & Peterson 



(|2006l) . which makes use of the FWHM of the Hfi emission-line and the continuum luminosity 
at 5100 A. Luminosity at 5100 A (L5100) vvas estimated from the median flux derived with 
the power-law fit done with STARLIGHT and the luminosity distances of the objects using 
the redshifts give n in NED. The bolometric luminosity was obtained using Lboi ~9/lL^5ioo (see 



Kaspi et al. 



20001), which yields the Eddington ratio Lboi/^Edd, where LEdd-l-SxlO^*^ (Mbh/Mq) 
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Our luminosity results are shown in Table |4l For objects without Hj 3 broad componen t, the Ha 



line FWHM was converted to Hfi FWHM using equation 3 given by 
mass estimates for our sample are shown in Table [51 



Shenetal 



(l2008h . The BH 



We can also derive the BH mass using the so called M-cTi, relation (|Ferrarese & Merritt 



20001 : 



Gebhardt et al. 



2000|) and see if these estimates are in agreement with the black hole mass 



derived for these objects using scaling relations. This is the reason why we obtained the velocity 



dispersion (ctj^) with STARLIGHT. Howeve r since the use of this relation depends on the bulge 

J.g.HE 



type associated with the host galaxy (e.g. lHull2008r) . we will present these estimates in Paper II, 
along with a detailed discussion of the bulge type present in each galaxy based on the photometric 
and spectroscopic observations. 



In Fig. [13] we show the BH mass vs. luminosity relation obtained for the objects in our 
sample. The diagonal lines show the Eddington limits shown by our objects. Note that internal 
extinction has not been taken into account. All objects have low to moderate accretion Eddington 
rates (Lboi/Ledd <0.11). 



4.2. Double-Peaked Objects 

Spectral analysis showed that among our objects we have five objects that clearly show 
narrow double-peaked emission lines: galaxies #1, #3, #4, #5 and #8 (c.f. Table [T])- In Table [3] we 
present the velocity differences, AV = Ked - V^biue, for Ha, [N II]/16548, [N II]/16583, [S U]A67\6, 
[S II]/16731, H;8,[0 III]i4959, and [O III]i5007 lines. Our definition of a narrow line 
double-peaked object requires that in all these lines a double-peaked Gaussian components should 
be present. The average values and standard deviation obtained for each source are listed in the 
bottom line of this table, which shows that they are within 1 14 < AV < 256 km s^^ 
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4.3. Notes on particular objects 

4.3.1. J120655.63+501737.1 (#1) 

The profile decomposition shows that this AGN was best fitted using double-peaked lines for 
its narrow line region (NLR, see Fig.|7]). Its spectrum was fitted with a broad component for the 
Hfi and Ha lines, and from its locus in the BPT diagrams shown in Fig. [121 we classify this object 
as a composite AGN. The blue component indicates that it is a Sy 1.8-l-SB, and the red component 
that it is a Sy 1.5 + SB. 



4.3.2. J121600.04+124114.3 (#2) 



In 



Mazzarella & Balzanol ( 



Veron-Cettv & VeronI (|2010() and in NED this object is classified as a Sy 1 .9, and in 
1" 



1986|) it is classified as a Starburst (SB) galaxy. We agree with both 



classifications since this object appears in one diagram as a composite Sy 1.9 + SB, and in the 
other as a clear Starburst galaxy, c.f., Fig.[l2l We noted that the emission line fitting of this object 
clearly needed a double-peaked Gaussian component for [O III], [S II] and [N II] lines; but when 
fitting the narrow components of Ha and HJ3 only one Gaussian component was required. Since 
we did not needed double-peaked Gaussians to fit the narrow line components of Ha and Hfi, this 
object was not classified as double-peaked. 



4.3.3. J121607.08+504930.0(#3) 

The profile decomposition analysis shows that this object is a narrow-line double-peaked 
AGN, see Fig. [8l Its spectrum was fitted with a broad component for the HJ3 and Ha lines, and 
from the locus of each component in the BPT diagnostic diagrams shown in Fig. [T^l we classify 
this object as a Sy 1.8 galaxy. 
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4.3.4. J141238.14+391836.5 (#4) 

From our profile decomposition, we find that this AGN is a narrow-line double-peaked 
(Fig. |9l) source. In addition, we noted that this object lacks a broad component in the Hfi region. 
From its locus in the BPT diagrams in Fig. [TH and using the blue and the red components, we 
classify this object as a Sy 1 .9 galaxy. 



4.3.5. J143031. 18+524225.8 (#5) 

This object is also a narrow-line double-peaked source, see Fig.[10l Since we need a broad 
Gaussian component to fit the Hor region, and from the BPT diagnostic diagrams in Fig.[T2l using 
both the blue and red components we classify this object as a Sy 1.5-1- SB. 



4.3.6. J 162952.88+242638.3 (#8) 

The [O III]/l/l4949,5007 doublets clearly show a double-peaked profile and a rather broad 
wing. We fitted the profile with two sets of three Gaussians. Two Gaussian profiles were used for 
fitting the narrow component of HyS (Fig.fTTI). The profile decomposition shows abroad Gaussian 
component for Ha. The BPT diagnostic diagrams in Fig. [T2l clearly show that, using both red 
and blue components, the object is a composite AGN and we classify it as a Sy 1.9-1- SB. The 
results on the spectral analysis of the XMM-Newton data on this source (see Appendix [A]) are 
in agreement with the optical data ones. This is based on the obtained value of the F index and 
the low values found for the equivalent Nr. Both values are compatible with Type I objects (e.g.. 



Piconcelli et al. 



2005 ) 
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5. Discussion and Conclusions 

The sample presented in this work consists of ten nearby intermediate-type AGN, some of 
them showing a weak broad Ha line component. An example of this was presented in the profile 
analysis of J162952. 88+242638. 3. It clearly shows a rather weak broad component in Ha (cf. 
Fig.fTTI). typical of a Sy 1.9 galaxy. This result is also supported with our analysis of the XMM 
public available data on this object that showed that the power index and the hydrogen column 
absorption are compatible with Type 1 (cf., Appe ndix A) objects. Th us, our Sy 1.9 classification 



for this object resulted to be in disagreement with 



Trippe et al. 



(|2010[) . that recently classify this 



object as a Seyfert 2 using the same optical and XMM data. 

The BPT diagnostic diagram [O III]/H;S vs. [S II]/Ha clearly shows that all single-peaked 
objects are composite AGN. In the case of the narrow line double-peaked objects detected in 
this work, three are composite AGN (#1, #5, and #8) and only two of them (#3 and #4) are pure 



intermediate Sy galaxies. Therefore, 80% 



+7.2% 



■17 



AGN. It is interesting to note that using the 



of our inter mediate-type objects are composite 



Whittle! (|l992h classification for intermediate Sy 



galaxies, we confirm that all objects studied in this work are intermediate-type AGN. 

The FWHM of HjS component was used to estimate the BH mass in all objects for 
the first time. We found that the sample studied in this work has a BH mass range of 
6.54+0.1 6 < log Mrh < 7.8 1 +0.14 . In particular three objects (#5, #9, and #10) have Mbh -IO^Mq. 



Recently, 



Kormendv et al. 



(120111) suggested that small black holes in bulgeless galaxies or 
pseudobulges grow as low-level Sy. This idea will be explored in our companion Paper II, where 
we study the bulge properties of this sample. Also, in Paper II the Mbh estimates obtained in 
this work will be compared with estimates for the same sample through the use of the M-cr+ 
correlation. We also found that objects in our sample have low to moderate accretion rates that go 
from (0.21+0.23)xl0"2 to (3.44+0.24)xl0"^ and are located at the low-tail of bolometric AGN 



luminosities. 
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Among the intermediate-type AGN identified in this sample, 50 %+ 14.4% show n a rrow 



Comerford et al. 



(l2009h 



double-peaked emission lines. From the BPT diagrams and following 
we find that two objects #3 and #4 are good candidates for dual AGN since both the blue and 
red emission line components are dominated exclusively by AGN activity. The other three 
narrow-double peaked objects #1, #5 and #8 show one emission line component in the AGN 
region while the other appears in the composite AGN -I- SB region. Therefore, the precise 
classification of each component in the case of double-peaked AGN is essential, since it gives 
hints on the origin of the observed narrow double-peaked lines. This information could be useful 
in identifying new dual AGN. 

The objects showing narrow double-peaked emission lines found in this work need to 
be studied through observations enabling spatially resolved imaging and spatially resolved 
spectroscopy. In searches of dual AGN, Type 1, Type 2 and intermediate-type AGN should be 
carefully separated, due to the high serendipitously number of narrow double-peaked sources 
found in our sample. These objects could be useful for establishing the properties and origin 
of narrow double-peaked emission lines profile in cases where the peak velocity separation lies 
below 300kms~^ 
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A. XMM DATA ANALYSIS OF J162952.88+242638.3 (#8) 

We searched for X-ray data available for our sample objects in Chandra and XMM-Newton 
databases. We found that only one of our objects has been observed with XMM-Newton. This 
object is 1162952.88-^242638.3, also known as Mrk 883. The three XMM-Newton observations 
analyzed for this object were performed during August 13th, 15th and 21st 2006 (Obs. ID 
0302260101; 0302260701 and 0302261001, respectively). The exposure times after flare 
removal (see below) were 7.9, 10.5 and 10.5 ks. AH the observations w ere processed using 



the standard Science Analysis System, SAS vlO.3.0 (iGabriel et al. 



20041) and using the most 



updated calibration files available in January 201 1. The EPIC-pn events were filtered due to high 



background events using the method described in 



Piconcelli et al. 



(|2004|) . No signs of pile-up or 



out-of-time events were detected in any of the observations. 

The spectra of each observation were extracted from circles centered in the maximum 
emission pixel and with radii which maximized the signal-to-noise, i.e., 22.5 ". The background 
circular regions used for the spectral analysis were located close to the source and free of any 
contaminating source. Both response matrices were generated for each observation using arfgen 
and remgen SAS ( [http : //xmm . esac . esa . int/ sas/[ ) tasks. The spectr a were binned to have a t 



1977h 



least 20 counts per bin in order to be able to use the modified technique (iKendall & Stuart 
for the data fit. All the three spectra were analyzed simultaneously. The shape of the continuum is 
quite similar among the three observations, with a slight increment of the emission for the third 
observation. After taking into account the absorption due to our Galaxy, Nh=3.S3 x lO^^cm"^, the 
continuum emission can be satisfactorily fitted with a single absorbed power-law with a photon 
index r=1.7+0.03 and an equivalent hydrogen column density of the order of 1 x lO^^cm"^. 
Table [6] summarizes the values of the parameters for the best fit model. 

None of the observations of this object present an iron line statistically significant, although 
both a neutral iron line with the energy fixed to 6.4 keV and an iron line with the energy left 
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free to vary were tested. During the eight days in which the observations span, moderate flux 
variability between the two first observations and the last one was detected. Neither power-law 
index nor intrinsic absorption variation were measured. Figure [T?l summarizes the results for the 
X-ray analysis of SDSS J162953. 88-^242638. 3: the best-fit model for the three spectra with the 
residuals to the data, and the contour plots for the confidence level between the Y index and the 
equivalent Nh for the three observations. 
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Fig. 1.— STARLIGHT fit of the SDSS spectra. The black line shows the original SDSS de- 
redshifted spectrum. The red line shows the pure stellar continuum spectrum and the purple line 
the power-law used to fit the AGN continuum. The blue line shows the best fit to the data and 
represents the sum of the red spectrum plus the AGN continuum. We show in green the final 
pure AGN spectra of J120655.63+501737.1 and J121600.04+1241143. Bottom panels show the 
normalized residuals. 
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Fig. 2.— Spectra of J121607.08+504930.0 and J141238. 14+391836.5. As in Fig.[IJ 
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Fig. 3.— Spectra of J143031. 18+524225. 8 and J144049.35+505009.2. As in Fig.m 
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Fig. 4.— Spectra of J153810.05+573613.1 and J162952.88+242638.3. As in Fig.[IJ 
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Fig. 5.— Spectra of J212851. 19-010412.4 and J234428. 81 + 134946.0. As in Fig.m 
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Fig. 6. — A blow-up of our best model obtained from Monte Carlo simulations for the the spectrum 
of J212851. 19-010412.4. This fit is shown also in blue in the left panel of Fig.[5l The lines used to 
estimate the velocity dispersion are marked with vertical dotted lines. 
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Fig. 7. — Spectral decomposition for SDSS J120655. 63+501737.1. Left and right panels show our 
best fit obtained for the blue and red part of the spectrum, respectively. In both panels black lines 
correspond to the pure AGN spectrum obtained with STARLIGHT, blue lines show the Gaussian 
components used for the fit and red lines the best fit obtained. At the bottom of each panel the 
corresponding normalized residuals are shown. The [O III]/14959,5007 lines and the narrow com- 
ponent of H/3 were best fitted using 2 Gaussian components. For the broad component of H/3 only 
one Gaussian was needed for the fit. The red part of the spectrum was modeled in a similar way, 
i.e., two Gaussian components for the narrow lines and one Gaussian for the broad component of 
Ha. 
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Fig. 8.— Spectral decomposition for SDSS J121607.08+504930.0. We have modeled the lines in 
a similar way as described in Fig. Ul 
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Fig. 9.— Spectral decomposition for SDSS J141238. 14+391836.5. This object has also been mod- 
eled using two Gaussian components for the narrow lines. But, when fitting the broad components 
we found no evidence for a broad component in Hfi. For the red part of the spectrum a broad Gaus- 
sian component was needed for properly fitting the Hor line, and the two Gaussian components for 
the narrow lines. 
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Fig. 10.— Spectral decomposition for SDSS J143031. 18+524225.8. This object has been modeled 
similar to the one shown in Fig. |9l 
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Fig. 11.— Spectral decomposition for SDSS J162952. 88+242638.3. The lines were fitted using 
three Gaussian components for [OIII]/I4959,5007 since it shows a wing component in addition to 
the two narrow components. The red part of the spectrum was fitted using two Gaussian compo- 
nents for the narrow lines. There is no broad component in Hfi, but a broad component was fitted 
for Ha. 
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Fig. 12. — Location of the o bjects of the sa mple in two BPT diagnostic diagrams. Left side, the dot 



ted lines come from: (IKewley et al 



2003L K03). Some objects 



1 1200 IL KOI) and (iKauffmann et al. 
are composite AGN, i.e., they lie between KOI and K03 lines, meaning that they have a Starburst 
galaxy (SB) plus an AGN. Filled black dots mark the locus of single-peaked objects. Objects with 
narrow double-peaked emission lines are shown twice since we used the ratios obtained from their 
red (open squares) and blue (open triangles) components. Rig ht side, BPT d i agrarn showing that 



all objects are located above the LINER region accordingly with 
classification is presented in Table [51 
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Fig. 14.— Left column: EPIC-pn observed spectra of SDSS J162952.88+242638.3 (Mrk 883) 
together with the best-fit model and the residuals. More details of the model can be found in 
Table [6l Right column: contour plots of the index of the power law and the equivalent hydrogen 
column density for the three XMM-Newton observations of this object. Each line corresponds to 
a different observation date. 
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Table 1. Intermediate-type AGN Sample 



Galaxy 


SDSS 


Other 


z 






# 


ID 


Name 




(Mpc) 




1 


J120655.63+501737.1 


SBS 1204+5()5B 


0.062 


278.01 


16.20 


2 


J121600.04+ 1241 14.3 


Mrk764 


0.066 


296.78 


14.87 


3 


J121607.08+504930.0 


Mrk 1469 


0.031 


135.92 


15.27 


4 


J141238.14+391836.5 


NGC 5515 


0.026 


113.58 


13.78 


5 


J143031. 18+524225.8 


SBS 1428+529 


0.044 


194.77 


15.25 


6 


J144049.35+505009.2 


SBS 1439+510 


0.117 


544.65 


17.28 


7 


J153810.05+573613.1 


SBS 1537+577 


0.073 


329.87 


15.51 


8 


J162952.88+242638.3 


Mrk 883 


0.038 


167.47 


15.26 


9 


J212851. 19-010412.4 


IC 1385 


0.052 


231.51 


15.21 


10 


J234428.81+134946.0 




0.068 


306.21 


16.32 



Note. — Col. 1: galaxy identification. Col. 2: SDSS id. Col. 3: other names. Col. 4: 
redshift by NED. Col. 5: luminosity distance. Col. 6: Petrosiang-band magnitude given 
in SDSS-DR7. 



Table 2. Optical Spectroscopic Parameters 



Galaxy S/N FWHM (Hff) FWHM (H/?) log[OIII]/H8 log[NII]/HQr log[SII]/Ha 

# (kms"') (kms"') blue red single-peak blue red single-peak blue red single-peak 



1 


13.5±0.4 


2950±58 


2565±284 


0.96±0.06 


0.20±0.01 




-0.31±0.02 


-0.17±0.10 




-0.54±0.03 


-0.57±0.03 




2 


9.2±0.2 


4860±158 


5549±449* 






0.18±0.01 






-0.34±0.02 






-1.14±0.05 


3 


8.8±0.2 


6700±69 


6190±739 


1.06±0.04 


1.24±0.05 




-0.13±0.01 


0.19±0.01 




-0.41 ±0.02 


0.23±0.01 




4 


8.3±0.2 


3551 ±279 


4068±489* 


0.35±0.01 


0.34±0.01 




-0.06 ±0.002 


-0.07±0.002 




-0.34±0.01 


-0.38±0.01 




5 


5.9±0.1 


3878±341 


4438±583* 


0.12±0.01 


0.37±0.02 




-0.19±0.01 


-0.13±0.01 




-0.53±0.02 


-0.52±0.02 




6 


6.5±0.4 


3194±143 


2584±519 






0.32±0.04 






-0.15±0.02 






-0.51 ±0.06 


7 


10.6±0.3 


2491±38 


2788±187 






0.16±0.01 






-0.34±0.02 






-0.64±0.04 


8 


11.5±0.3 


2776±87 


3187±219* 


0.76±0.04 


0.18±0.01 




-0.29±0.02 


-0.35±0.02 




-0.12±0.01 


-0.48±0.02 




9 


10.7±0.4 


1955±25 


2443±101 






0.42±0.03 






-0.24±0.02 






-0.66±0.05 


10 


6.3±0.5 


2094±38 


2442±178 






0.52±0.10 






-0.36±0.07 






-0.70±0.13 



Note. — Col. 1: Galaxy # (cf. Table[T). Col. 2: SPSS spectra S /N obtained from STARLIGHT Cid Femandes etauLoOS ). Col. 3: FWHM(Hq-) broad component. Col. 4: FWHM(HyS) 
broad component. Cols. 5-7: log [O III]/H/3 line ratios for double-peaked objects (blue and red) and for single-peaked sources. Col. 8-10:log [Nin]/H^ line ratios for double-peaked objects (blue 
and red) and for single-peaked sources. Col 11-13: log [SII]/Hq' line ratios for double and single-peaked sources. Note that we mark an * in objects where the FWHM(H^) was obtained using 



Shen et al 



2008 ), see text. 



Table 3. Velocity differences in narrow double-peaked AGN 



SDSS 


J120655.63+501737.1 


J121607.08+504930.0 


J141238.14+391836.5 


J143031. 18+524225.8 


J162952.88+242638.3 




(#1) 


(#3) 


(#4) 


(#5) 


(#8) 




AV 


^v 


AV 


AV 


AV 


Line 


(kms~') 


(km s"') 


(kms"') 


(km s'') 


(km s"') 


HQf^6562.74 


116±17 


209±122 


264±15 


296±5 


115±3 


[N 111^6548 


117±31 


177±48 


232±28 


276±32 


110±17 


[NII]i6583 


118±5 


160±107 


231 ±32 


274±12 


135±5 


[S II]i6716 


111±15 


97±9 


226±7 


272±15 


128±7 


[SII]/16731 


96±35 


56±8 


226±13 


271 ±20 


123±7 


1^3^4861.29 


120±11 


160±10 


263±80 


264±17 


142±20 


[Oin]a4959 


118±35 


132±26 


258±18 


199±185 


139±10 


[Oin]a5007 


116±11 


120±7 


256±120 


198±64 


130±3 


average value 


114±7 


139±45 


245±16 


256±34 


128±11 



Note. — Velocity differences obtained between the red and the blue Gaussian components (AV = Vred - Vbiue)- Col. 1: line rest frame. Cols. 
2-6: velocity differences obtained for SDSS objects J120655.63+501737.1, J121607.08+504930.0, J141238.14+391836.5, J14303 1.1 8+524225.8, 
and J162952.88+242638.3, respectively. Last row indicates the average AV value and standard deviation of each object considering all lines 



Table 4. Luminosities and Eddington Rates 



Galaxy 




Lbol 




Lbol/LEdd 


# 


(10*" erg s"') 


(10''2 ergs-l) 


(lO'^ergs-l) 


(10-2) 


1 


49.92±6.98 


45±6.28 


15.05±5.08 


2.98±1.42 


2 


14.52±4.84 


13.07±4.36 


37.99±14.24 


3.44±0.24 


3 


42.02±3.65 


37.82±3.28 


80.45±26.41 


0.47±0.20 


4 


5.77±1.53 


5.19±1.38 


12.87±5.39 


0.40±0.28 


5 


2.11±0.41 


1.90±0.37 


9.27±5.38 


0.2I±0.16 


6 


40.43±6.83 


36.39±6.14 


13.77±7.32 


2.64±0.02 


7 


45.47±4.90 


40.92±4.37 


16.98±3.96 


2.41±0.82 


8 


53.74±7.17 


48.37±6.45 


24.12±6.03 


2.01±0.01 


9 


18.14±1.63 


16.33±1.47 


8.23±1.43 


1.98±0.52 


10 


5.08±1.76 


4.57±1.59 


4.35±1.59 


1.11±0.75 



Note. — Col. 1: galaxy # (cf. Table [T). Col. 2: luminosity at 5100 A 
measured using a continuum power-law fit. Col. 3: bolometric luminosities 
obtained assuming Li,oi = 9/iL^5ioo- Col. 4: Eddington luminosities using 
our estimated MgH. Col. 5: Eddington ratio. 



Table 5. BH Mass 



Galaxy 




log Mbh 


Spectral Classification 


Double-peak 


# 


(km s-') 


(Mo) 






1 


90±6 


7.08± 0.14 


Sy 1.8 + SB/Sy 1.5 + SB 


yes 


2 


111±6 


7.48±0.16 


Sy 1.9 + SB/SB 


no 


3 


147±7 


7.81±0.14 


Sy 1.8 


yes* 


4 


195±8 


7.01± 0.18 


Sy 1.9 


yes* 


5 


165±9 


6.86±0.25 


Sy 1.5 + SB 


yes 


6 


188±17 


7.04±0.23 


Sy 1.5 + SB 


no 


7 


91±6 


7.13±0.10 


Sy 1.5 + SB 


no 


8 


187±9 


7.28±0.11 


Sy 1.9 + SB 


yes 


9 


79±6 


6.82±0.08 


Sy 1.5 + SB 


no 


10 


93±11 


6.54±0.16 


Sy 1.8 + SB 


no 



Note. — Col. 1: galaxy id. (cf. Table[T). Col. 2: velocity dispersion obtained with 
STARLIGHT. Errors were estimated from our Monte Carlo simulation on the data, 
see text. Col. 3: Black hole mass estimates were obtained with the relation given by 



Vestergaard & Petersor 



200d) . Black hole mass for object #8 have to be taken with 
care since the host galaxy is probably a merger. Col. 4: Sy type. Col. 5: yes if it is 
double peak, and * shows the two dual AGN candidates found in our study. 



Table 6. X-ray Spectroscopic Parameters 



Obs. 




Neutral Abs 


Power Liiv^" 


Emission Line 




CiOodncss 




Flux 


Luminosity 


Day 


Model 


Nh 


r 


Energy 


EW 


xl 


d.o.f. 


F-test 


0.5-2 keV 


2-10 keV 


0.5-2 keV 


2-10 keV 


13 


A 
B 


9±2 
9±2 


1.68±0.06 
1.69±0.06 


6.4f 


<170 


134.0 
133.9 


118 
117 


35.2% 


65+°'" 

"•"■^-0.04 


1.51 ±0.05 


2 64+012 


4.90 ±0.16 


15 


C 
A 


10±2 

9±2 


1.70±0.06 

1.71±0.05 


6.88+0;?^ 


210+180 

^^"-160 


129.0 
135.0 


116 

129 


87.5% 


0.66±0.02 


1.47±0.05 


2.66±0.08 


4.79±0.16 




B 


10±2 


1.72±0.05 


6.4f 


210+'° 

^^"-150 


128.9 


128 


98.5% 








U) 


21 


C 
A 


10±2 
11±2 


1.72±0.05 
1.69±0.05 


-0.10 


200+13" 


127.3 
128.9 


127 
145 


97.6% 


0.83±0.03 


1.95±0.05 


3.43±0.12 


1 

6.33±0.16 




B 


11±2 


1.70±0.05 


6.4f 


75+^ 


126.9 


144 


86.6% 












C 


11+2 


1.7()±().()T 


6.5^+0.15 


150+1 10 


123.9 


143 


94. 1 











Note. — Col 1: day of observations, data are from August 2006. Col 2: model A: absorbed power law; model B: absorbed power law plus a neutral iron line with 
the energy fixed; model C: absorbed power law plus an iron line with the energy left free to vary. Col 3: neutral absorption in units of lO^" cm"^. Co (4): power law 
index. Col (5) and (6): energy in keV and equivalent width in eV, respectively. Col (7): modified reduced Col (8) degrees of freedom. Col (9) F-test. Col. 10 and 
11: fluxes in lO^i^ erg s^i cm^^ Col. 12 and 13: luminosities in Vf^^x% s^i . 



